May, 1993 



IFUP-TH 17/93 



The Three-Loop Lattice Free Energy 

B. Alles a , M. Campostrini 6 , A. Feo & , H. Panagopoulos & ' c 

a ) Departamento de Ffsica Teorica y del Cosmos, Universidad de Granada, Spain 
b ) I. N.F.N, and Dipartimento di Fisica dell'Universita, Pisa, Italy 
c ) Dept. of Natural Sciences, Univ. of Cyprus, Nicosia, Cyprus 



ABSTRACT 

We calculate the free energy of SU(iV) gauge theories on the lattice, to 
three loops. Our result, combined with Monte Carlo data for the average 
plaquette, gives a more precise estimate of the gluonic condensate. 
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In a previous work we described a scheme for carrying out perturbative calculations 
on the lattice, using a symbolic manipulation language. We also used this scheme to obtain 
the first three-loop result, related to the topological susceptibility. 

In this paper, we present a three-loop calculation of the free energy in pure SU(iV) gauge 
theories. Using this perturbative result together with Monte Carlo data for the average pla- 
quette we arrive at a more precise determination of the (non- perturbative) gluon condensate 

This condensate has been widely used in phenomenological studies of QCD, in particular in 
the context of the SVZ sum rules @. In addition, the feasibility of extracting G<i from a 
lattice simulation is a challenging theoretical issue on its own right; the presence of mixing 
with lower dimensional operators {F^ V F^ V mixes with the unit operator) represents a priori 
a major stumbling block. It is important to see just to what accuracy this extraction can be 
pushed given the systematic error. 

The gluonic condensate is related to Monte Carlo data through the expression 

(n> MC = Z(p) {-'Qj G 2 a 4 + (n> pert , (2) 



where 



and IT is the plaquette operator 

The quantities Z((3) and (n) per t can be calculated perturbatively; while Z{j3) for the plaquette 
operator is very close to 1 and can be left out to a rather good approximation, (n) per t is far 
from negligible, and it is important to determine it with the best possible precision. 
We write: 

<nW = iv + ! + £| + J + ..., (4) 

where c\ = —N(N 2 — l)/4 and C2 was calculated in [Ref. 3| . 

We have carried out the calculation of the three-loop coefficient C3. This a priori entails 
computing all vacuum diagrams with one insertion of the plaquette operator; there are 63 
such distinct diagrams. The relationship between IT and the Wilson action 

S=J2( N ~^) (5) 

X,[IV 

allows a simplification of this task, by relating (n) per t to the partition function 
Z = f[dU] exp(-(3/2N S). We have: 



(n) P ert = I J [dU] e~^ S U x ,^ 



N- — - f[dU]e-^ s S = N + — ^-\nZ . 
12V Z J 1 J 12V d(3 



(6) 
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Now, In Z has the form 



which implies that C3 = —N/2d,2- The calculation of (I2 involves computing all 3- loop con- 
nected vacuum diagrams with no operator insertions; there are 29 diagrams giving nonzero 
contributions altogether, shown in Fig. 2. 

We have computed these diagrams using our algebraic package described in [Ref. 1| . The 
main building blocks of this package are: 

i) Generating all relevant vertices. 

ii) Performing contractions for each diagram and producing the corresponding integrand. 
Here, as well as in i), it is crucial to exploit all available symmetries, in order to render 
the resulting expressions more compact. 

iii) Integrating over loop momenta, for infinite as well as finite lattices. 

The contributions from each diagram in the Feynman gauge are given for reference in Table 
I. The final result, given in Table II, is, of course, gauge independent. Incidentally, the CPU 
time necessary for the actual runs (once the development and testing stages are over) is ~ 3 
hours on a Sun SPARC ELC to produce the integrand for some of the more complicated 
diagrams and ~ 2 hours for a highly optimized integration on lattices up to 16 4 . (A notable 
exception are diagrams with the topology of Fig. 21, known as "Mercedes-Benz emblem" ; the 
integration of these diagrams requires ~ 20 hours already for an 8 4 lattice.) 

A preliminary presentation of the results reported here was given in [Ref. 4j . Note however 
a discrepancy due to an error in the preliminary version. 

Using numerical results for L 4 lattices (L = 2, . . . , 16) we can extrapolate out to L = 00 
with a functional form: 

r(L) = r + £ (+^) i,j = 2,3,4 (8) 

depending on the diagram, and compare with the infinite lattice results r(oo). (Infinite 
lattice results are obtained with the Schwinger representation technique, explained in [Ref. 1| .) 
The discrepancy is typically a fraction of a per mille, so that r(L) can be used also for 
lattices with L > 16. The actual analytic form of r(L), for L — > 00, could also involve 
terms proportional to (In L) n /L m ||; however the introduction of these terms does not 
alter our extrapolations. It should also be mentioned here that, while lattice integrals are 
ultraviolet convergent, individual diagrams may exhibit infrared divergences, which cancel 
only when they are summed up. For the case at hand, the diagrams in Fig. 2o, containing a 
double insertion of the one-loop renormalized gluon propagator (Fig. 1), must be evaluated 
simultaneously, in order to guarantee infrared convergence. 

Having obtained the value of C3 we can now reanalyze Monte Carlo data for SU(2) 0, 
and SU(3) [§, using Eq. (2), to extract the gluon condensate. We perform two series of 
X 2 fits, one in which C3 is a parameter to be fitted, the other using the calculated value of C3 
and fitting instead C4. The results are listed in Table III; in all cases we have x 2 /d.o.f. ~ 1. 
It is reassuring to see that, despite variations in the fitted values of the parameters Cj, the 
condensate exhibits a very stable behavior. 

To conclude, our algebraic scheme has made possible this 3-loop calculation; done by 
hand, such a calculation would be extremely cumbersome, if not downright impossible. At 



3 



the present stage of development of this scheme, one can compute vacuum expectation values 
of any operator to an arbitrary number of loops (modulo computer limitations). We are 
currently extending the scheme to include matrix elements with external fields; the main 
problem there is extracting the analytic dependence on external momenta, which can be a 
rather delicate issue beyond one loop. We will report on this subject in a future publication. 
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Table I 

The contribution to C3 of individual diagrams, shown in Fig. 2, in the Feynman gauge. We 
use an L 4 lattice and gauge group SU(AT). Each entry must be multiplied by N 5 (N 2 — 1) 10~ 4 . 



Fig. 


L = 8 


L = 12 


L = 16 


L = 00 


a 


41.874 -93.004/iV 2 


41.833 - 92.971/7V 2 


41.811 - 92.946/7V 2 


41.775 - 92.902/7V 2 




+52.045/7V 4 


+52.076/7V 4 


+52.081/7V 4 


+52.083/iV 4 


b 


1.0803 


1.0974 


1.1035 


1.1113 


c 


38.582- 43. 395/iV 2 


40.013 - 44.714/7V 2 


40.528 - 45. 184/7V 2 


41.200- 45. 795/7V 2 


d 


1.0806 


1.0975 


1.1035 


1.1113 


e 


1.0803 


1.0974 


1.1035 


1.1113 


f 


0.22019 


0.23398 


0.23894 


0.24543 


g 


0.018161 


0.021014 


0.022139 


0.023730 


h 


-10.464 


-10.696 


-10.781 


-10.892 




+ {l/N 2 -3/N 4 ) 26.993 


+ (l/7V 2 -3/7V 4 ) 27.013 


+ (l/iV 2 -3/7V 4 ) 27.017 


+ (l/iV 2 -3/iV 4 ) 27.019 


i 


0.41123 


0.41825 


0.42074 


0.42394 


j 


-1.7753 


-1.7876 


-1.7880 


-1.7839 


k 


0.084525 


0.089997 


0.091975 


0.094566 


1 


-0.50731 


-0.56679 


-0.58988 


-0.62246 


m 


0.023925 


0.027060 


0.028245 


0.029881 


n 


0.036795 


0.037688 


0.038022 


0.038479 





-135.184+ 290.471/iV 2 


-136.433 + 291.877/7V 2 


-136.861 + 292.342/iV 2 


-137.404 + 292.919/iV 2 




-156.136/A^ 4 


-156.227/iV 4 


-156.243/iV 4 


-156.25/iV 4 


Total 


-63.438+ 181. 064/iV 2 


-63.517+ 181. 204/iV 2 


-63.532+ 181. 229/iV 2 


-63.537+ 181. 240/7V 2 




-185.068/7V 4 


-185.191/iV 4 


-185.211/iV 4 


-185.223/iV 4 



Table II 

The calculated, gauge independent value of C3 in SU(2) and SU(3). 



C3 


L = 8 


L = 12 


L = 16 


L = 00 


N=2 
N=3 


-0.28549 
-8.8655 


-0.28599 
-8.8781 


-0.28608 
-8.8805 


-0.28611 
-8.8814 



Table III 

We list the values of tt 2 /12N G 2 /Aq CD (for gauge groups SU(2) (0=2.45) and SU(3) 
(/3=6.0)), as obtained from Eq. (2) and Monte Carlo data through a series of fits, in which: 
a) C3 was an additional parameter to be fitted, or: b) the exact value of C3 was taken from 
our calculation, fitting instead c\. 



yr 2 /12iVG 2 /A 4 CD 


SU(2) 


SU(3) 


a 
b 


1.09(09) 10 7 
1.06(10) 10 7 


2.79(33) 10 8 
2.39(61) 10 8 
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Fig. 1 The one-loop gluon propagator. A square 
stands for the integration measure. 




